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半 透 明 介 质 场 参数 的 非 接触 测量 算法 研究 
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摘 要 : 本 文 在 考虑 半 透 明 介 oes 发 射 以 及 散射 衰减 作用 的 情况 下 采用 视 在 光线 法 计算 某 探测 方向 上 介质 边界 出 射 辐射 强 
度 分 布 作为 逆 问 题 的 输入 数据 ， | LSQR 方法 在 介质 辐射 特性 参数 已 知情 况 下 重建 介质 三 维 温度 分 布 ， 重建 结 果 显 示 ， 无 论 
有 无 测量 误差 存在 ， 利 用 LSQR Oe en 建 出 介质 的 三 维 温度 场 。 在 此 基础 上 提出 
LSQR-SPSO 混合 算法 ， 并 应 用 于 同时 重建 三 维 温度 场 及 辐射 特性 参数 (吸收 系数 、 散 射 系数 )。 计 算 结果 显示 ， 无 论 有 无 测量 
误差 存在 ， 利 用 LSQR-SPSO 混合 算法 都 可 以 很 好 地 同时 重建 出 介质 的 吸收 系数 、 散 射 系 数 及 三 维 温 度 场 ， 相 比较 而 言 ， 介 质 
的 温度 场 更 容易 被 重建 出 来 。 
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Abstract: For the semi-transparent media which only consider the absorption, emission and Scattering attenuation 
characteristics, the outgoing radiative intensities at the boundary surface were Simulated by the line-of-sight (LOS) 
method, which were served as input for the inverse analysis. The least-square QR decomposition (LSQR) algorithm 
was used to estimate the three dimensional temperature distributions of the semi-transparent media with known 
radiative properties. And the retrieval results show that the temperature distributions could be retrieved accurately 
even with noisy data. For the participating media with unknown radiative propetrties, a hybrid least-square QR 
decomposition (LSQR) - Stochastic Particle Swarm Optimization (LSQR-SPSO) algorithm was developed to 
estimate the 3-D temperature distributions and radiative properties (absorption coefficients and scattering 
coefficients) simultaneously. And the retrieval results show that the absorption coefficients, scattering coefficients 
and 3-D temperature distributions also could be retrieved accurately even with measurement errors. It was also found 
that the temperature field could be estimated more accurate than the radiative properties. 

Key words: LSQR-SPSO; Temperature field; Radiative properties; Inverse problem 


0 引言 


目前 对 于 介质 温度 的 测量 主要 分 为 接触 式 测量 。 就 是 建立 高 效 、 准 确 、 稳 定 的 介质 温度 场 及 物性 场 
和 非 接触 式 测量 两 类 。 接触 式 测量 主要 是 指使 用 各 重建 算法 ， 这 也 成 为 近年 来 国内 外 学 者 广泛 关注 的 


种 不 同类 型 热电 个 、 示 温 漆 、 热 电阻 等 接触 式 物 理 研究 热点 。 刘 林 华 、 谈 和 平等 0 提出 了 网 络 搜索 与 


探 针 进行 测量 ， 但 这 些 经 典 的 接触 式 测量 方法 不 仅 ee 、 
干扰 流 场 ， 测 量 误差 大 ， 而 且 在 很 多 场合 无 法 使 用 ， 轿 梯 度 相 结 合同 时 反 演 一 维 半 透 明 介质 内 辐射 源 项 


更 为 重要 的 是 接触 式 测量 方法 很 难 实现 温度 分 布 测 。。 和 壁面 黑 度 的 混合 方法 。 齐 宏和 阮 立 明 等 ”首次 将 
量 。 非 接触 式 辐射 成 像 测 温 技术 通过 对 介质 辐射 图 ”智能 微粒 群 引 入 辐射 源 项 和 辐射 物性 的 逆 问 题 研 究 
像 处 理 进行 温度 测量 ， 具 有 测量 精度 高 、 非 插入 式 、 。 中 ， 发 现 该 算法 具有 对 不 依赖 初 值 且 无 需求 解 目标 
实时 连续 测量 等 优点 ， 且 能 够 获得 三 维 温度 分 布 函数 导数 等 优点 。 周 怀 春 、 郑 楚 光 等 中 提出 了 联合 反 
息 ， 无 多 代表 着 一 种 非常 有 前 景 的 技术 研究 站 势 。 演 一 维 介质 温度 和 辐射 物性 的 方法 。 娄 春 等 加 提出 了 
但 该 技术 不 能 通 通过 CCD 相机 或 热 像 仪 测试 技术 直接 一 种 解 耦 方法 用 于 联合 重建 一 维 光 学 薄 火 焰 或 非 散 
获得 介质 内 部 温度 分 布 ， 需 要 建立 复杂 的 辐射 道 问 本 aed) 
le ee 射 性 弥散 介质 的 温度 分 布 和 辐射 物性 参数 。 刘 冬 等 
题 计 算 模型 ， 且 通常 需要 已 知 介质 的 辐射 物性 参数 。 。 一 本 
对 于 非 接触 式 辐射 成 像 测 温 技 术 ， 其 核心 问题 。 于 二 县 仅 一生 QR 分 解法 TSQR》 全 全 全 国生 
, ”” ”“ “ 洛 法 对 二 维 参 与 性 介质 温度 场 和 辐射 物性 参数 进行 
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反 演 研究 。Klose 等 中 基于 有 限 差分 离散 坐标 法 和 准 
牛顿 法 研究 了 稳 态 传输 光学 成 像 问题 。 综 上 所 述 ， 
目前 对 于 参 de 度 场 重建 研究 ， 多 为 一 维 、 
二 维 温度 场 的 检测 ， 多 数 需要 已 知 或 者 假定 介 
质 辐射 特性 参数 ， We 
同时 重建 研究 较 少 。 
与 性 介质 温度 场 重建 问题 是 典型 的 不 适 定 问 
题 ， Tikhonov 正则 化 方法 、 截 断 
奇异 值 分 解 (TSVD ) 方法 及 LSQR 方法 等 。 本 文 将 
针对 考虑 吸收 、 发 射 及 散射 衰减 的 半 透 明 介 质 ， 在 
辐射 特性 参数 已 知 条 件 下 采用 LSQR 方法 利用 探测 
器 在 介质 边界 处 接收 到 的 辐射 强度 信息 重建 介质 三 
维 温度 分 布 。 当 介质 辐射 特性 参数 未 知 时 ， 提 出 随 
机 微粒 群 (SPSO ) 算法 与 LSQR 算法 结合 的 
LSQR-SPSO 混合 算法 ， 对 半 透 明 介 质 三 维 温度 场 和 
辐射 特性 参数 吸收、 散射 系数 ) 进行 同时 重建 。 


1 辐射 传输 正 问题 

在 仅 考虑 半 透 明 介质 的 吸收 、 发 射 以 及 散射 误 
减 作用 时 可 采用 视 在 光线 法 计算 任 一 探测 方向 上 的 
辐射 强度 ， 此 时 的 辐射 传输 方程 可 表示 成 如 下 形式 : 


0 


ED ,jer) = (0) (0) (1) 
式 中 ，1(r,t) 为 方向 辐射 强度 ; z 为 光学 厚度 ， 
T=K:S; Kk, 为 吸收 系数 ; wo 为 反照 率 ， 


oO=A / (Ka + K, ); 


ks 为 散射 系数 。 


1 几何 模型 
Fig. 1 The geometric model 
考虑 如 图 1 所 示 的 圆柱 形 半 透 明 介质 模型 ， 四 
周 为 透明 边界 条 件 ， 忽 略 环境 辐射 , 则 针对 任意 方 
向 的 边界 出 射 辐射 强度 表达 式 为 ， 


1,(r)= 7 (1 —exp(—t, )) 
nl (2) 
和 I -on > jc di)7w 


j=i+l 

式 中 ， n 为 r 方向 某 一 探测 线 所 经 过 的 网 格 单元 数 
目 ; zj; 为 探测 线 经 过 第 j 个 网 格 单元 的 光学 厚度 ; 必 
为 探测 线 经 过 第 i 个 个 网 格 单元 的 自身 黑体 辐射 强度 ， 


下 =a7 :了 7 为 探测 线 经 过 第 ;个 网 格 单元 的 温度 。 

将 r 方向 的 探测 平面 划分 为 M 个 单元 ， 则 可 得 

到 介质 内 任意 单元 温度 与 介质 边界 出 射 辐射 强度 之 
间 的 方程 组 ， 写 成 矩阵 形式 如 下 ， 

AT = (3) 

式 中 ，A 为 系数 矩阵 ; 也 为 黑体 辐射 强度 向 量 ; 工 为 

出 射 辐射 强度 向 量 ， 即 探测 器 接收 到 辐射 强度 向 量 。 


2 温度 场 及 辐射 特性 参数 重建 方法 


2.1 基于 LSQR 方法 重建 三 维 温度 场 

在 辐射 特性 参数 已 知 条 件 下 ， 方 程 组 (3) 为 线性 
方程 组 ， 因 此 可 直接 采用 LSQR 方法 求解 式 (3)， 由 
探测 器 接收 到 的 辐射 强度 分 布 重建 介质 三 维 温 度 
场 。 LSQR 方法 是 Paige 和 Saunders 四 于 1982 年 提出 
的 一 种 特别 适合 求解 系数 矩阵 为 大 型 、 稀 玻 和 窍 阵线 
性 方程 组 的 方法 。LSQR 方法 求解 的 思路 是 首先 把 任 
意 系 数 和 矩阵 方程 化 为 系数 矩阵 为 方 阵 的 方程 ， 然 后 
利用 Lanczos 方法 ， 求 解 方程 的 最 小 二 乘 解 。 


2.2 基于 LSQR-SPSO 混合 算法 同时 重建 三 维 温度 场 
及 辐射 特性 参数 

若 介 质 辐射 特性 参数 未 知 ， 则 需 根据 探测 器 接 
收 到 的 辐射 强度 分 布 同 时 重建 介质 温度 场 及 辐射 特 
性 参数 ， 此 时 方程 组 3) 为 非 线性 方程 组 ， 不 能 采用 
LSQR 方法 直接 求解 。 因 此 ， 本 文 引入 SPSO 算法 ， 
并 将 其 与 LSQR 方法 结合 同时 重建 介质 的 三 维 温度 
场 及 辐射 特性 参数 。SPSO 算法 最 早 是 由 Zeng 等 外 
提出 的 能 够 保证 以 概率 1 收敛 于 全 局 最 优 解 的 算法 ， 
广泛 应 用 于 各 种 优化 问题 的 求解 中 。 与 基本 PSO 入 
法 相 比 ，SPSO 去 掉 了 微粒 先前 的 速度 项 ， 使 得 速度 
失去 记忆 性 ， 从 而 减弱 了 全 局 搜索 能 力 ， 但 这 样 也 
使 得 在 进化 的 每 一 代 均 至 少 有 一 个 微粒 由 于 处 于 微 
粒 群 的 历史 最 好 位 置 而 停止 进化 ， 利 用 停止 进化 的 
微粒 改善 全 局 搜索 能 力 是 SPSO 算法 的 基本 思想 。 

对 于 温度 场 和 辐射 特性 参数 同时 重建 问题 需要 
两 个 不 同 探测 方向 如 和 yw, 的 探测 器 接收 到 的 辐射 
强度 分 布 I 和 I，, 作为 已 知 条 件 ,利用 SPSO 算法 搜 
索 吸 收 系数 ， 采 用 LSQR 方法 根据 y 方向 的 辐射 强 
度 分 布 I 计算 假设 的 温度 场 T', 然后 根据 假设 的 温 
度 场 T' 由 正 问 题 计算 得 到 假设 的 y, 方向 的 辐射 强 
度 分 布 L, ,并 与 wy, 方向 探测 器 探测 到 的 真实 的 辐射 
强度 I，, 构造 目标 函数 ， 如 式 (4) 所 示 ， 利 用 SPSO 算 
法 搜索 吸收 系数 使 目标 函数 极 小 化 ， 进 而 同时 重建 
出 介质 的 吸收 系数 、 散 射 系数 以 及 三 维 温度 分 布 。 


工 ， 一 工 
| n2 n2| (4) 


| 
n 


3 重建 结果 及 分 析 

如 图 1 所 示 的 圆柱 形 几 何 模型 ， 半 径 为 
r=0.5m ， 轴 向 长 度 为 L=20m ， 网 格 划 分 为 
No xNrxNz=1x10x20， 探 测 线 数 为 M =10000 ， 
按 (5) 式 假设 介质 的 真实 温度 场 ， 


T(z,r)=- (2- 20.0)x(0.5—3.0r? +2.0r? )+800 (5) 
式 中 ，z 为 圆柱 体 轴 疝 坐标 ，7 为 圆柱 体 径 向 坐标 。 


3.1 基于 LSQR 方法 的 温度 场 重建 结果 分 析 

首先 假设 介质 的 辐射 特性 参数 是 已 知 的 ， 假 设 
吸收 系数 为 k, =0.8 ， 散 射 系数 为 k, =1.8 ， 介 质 真 
实 温度 分 布 如 式 (5) 所 示 ， 分 别 在 无 测量 误差 以 及 
1%、3% 测 量 误 差 的 条 件 下 采用 LSQR 方法 根据 探测 
方向 2=90"，g =45” 上 的 探测 器 接收 到 的 辐射 强度 
分 布 重建 介质 三 维 温 度 分 布 ， 重 建 结 果 如 图 2 所 示 。 


Z/m 
(c) 7 =3% 


图 2 辐射 特性 参数 已 知 时 温度 场 重建 结果 


Fig.2 The reconstruction results of temperature field with known 
radiative properties 

从 图 2 中 可 以 看 出 ， 无 论 有 无 测量 误差 都 可 以 

得 到 很 好 地 温度 场 重建 结果 ， 不 同 测量 误差 条 件 下 

的 各 网 格 单元 重建 误差 分 布 如 图 3 所 示 ， 从 图 中 可 

以 看 出 在 无 测量 误差 时 温度 场 的 重建 误差 分 布 在 
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Relative retrieval error of temperature 
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图 3 辐射 特性 参数 已 知 时 温度 场 重 建 误差 分 布 


Fig.3 The reconstruction error distribution of temperature field 


with known radiative properties 


3.2 辐射 特性 参数 和 温度 场 同 时 重建 结果 分 析 

假设 介质 的 辐射 特性 参数 未 知 ， 设 吸收 系数 真 
值 为 x, =0.4m!， 搜 索 范 围 为 [0.1m1~1.0m!]， 散 射 
系数 真 值 为 x, =1.6m', 搜索 范围 为 [1.0m 1~2.0m1]， 
介质 真实 温度 分 布 如 式 (5) 所 示 ， 分 别 在 假设 无 测量 
误差 以 及 测量 误差 分 别 为 1%、3% 的 条 件 下 采用 
LSQR-SPSO 方法 根据 探测 方向 为 60=90*，@g=35? 
和 0=60°*，gp=45? 两 个 方向 上 的 探测 器 接收 到 的 辐 
射 强度 分 布 同时 重建 介质 吸收 系数 、 散 射 系 数 及 三 
维 温 度 分 布 。 不 同 测量 误差 条 件 下 的 吸收 系数 和 散 
射 系 数 重建 结果 如 表 1 所 示 。 

对 于 辐射 特性 参数 的 反 演 结果 ， 从 表 1 中 可 以 
看 出 ， 在 无 测量 误差 的 情况 下 ， 介 质 的 吸收 系数 和 
散射 系数 可 以 被 精确 地 反 演 出 来 ， 而 存在 测量 误差 
时 介质 吸收 系数 和 散射 系数 的 重建 误差 随 测量 误差 
的 增 大 而 增 大 ， 但 都 在 合理 的 范围 内 。 因 此 无 论 是 
否 存在 测量 误差 , 采用 LSQR-SPSO 混合 算法 都 可 以 
将 介质 的 吸收 系数 和 散射 系数 合理 地 重建 出 来 。 

对 于 温度 场 的 反 演 结果 ， 不 同 测量 误差 条 件 下 
A eel ei 
以 看 出 在 无 测 车 误差 时 温度 场 的 重建 误差 分 布 在 
10” 量 级 ， 随 着 测量 误差 的 增 大 ， 重 建 误差 也 随 之 增 
大 ， 即 使 存在 3% 的 测量 误差 的 时 其 温度 场 重建 误差 


107~102 量 级 ， 随 着 测量 误差 的 增 大 ， 重 建 误 差 也 
随 之 增 大 ; 但 即使 存在 3% 的 测量 误差 时 其 温度 场 重 
建 误 差 最 大 值 只 有 0.96%， 其 重建 误差 平均 值 为 
0.12%， 因 此 说 明 在 辐射 特性 参数 已 知 的 情况 下 ， 使 
用 LSQR 方法 可 以 很 好 地 重建 出 介质 的 三 维 温度 场 ， 
该 方法 具有 较 好 的 计算 精度 和 计算 稳定 性 。 


最 大 值 只 有 1.26%， 其 重建 误差 平均 值 为 0.70%。 由 
此 可 以 看 出 采用 LSQR-SPSO 算法 可 以 较 好 地 同时 
重建 出 介质 的 辐射 特性 参数 和 三 维 温度 场 。 

同时 发 现在 对 辐射 特性 参数 和 温度 场 同 时 重建 
时 ， 温 度 场 很 容易 被 重建 出 来 ， 并 具有 较 高 的 精度 ， 
而 辐射 特性 参数 重建 结果 相对 较 差 ， 较 难 重 建 出 来 ， 
这 可 能 是 由 于 测量 数据 对 辐射 特性 参数 的 敏感 性 较 
小 或 者 是 辐射 特性 参数 存在 多 值 性 造成 的 。 


| 


表 1 介质 吸收 系数 和 散射 系数 重建 结果 


Table 1 The reconstruction results of absorption coefficient and scattering coefficients 


y=0 
Inverse parameter True value 
Retrieval value  é&.% 


x, [oo] 0.4 0.4001 0.0250 
K [oo] 1.6 1.5999 0.0063 
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图 4 辐射 特性 参数 未 知 时 温度 场 重 建 误差 分 布 


Fig. 4 The reconstruction error distribution of temperature field 


with unknown radiative properties 


4 结论 

本 文 在 仅 考 虑 半 透 明 介质 吸收 、 发 射 以 及 散射 
衰减 作用 的 情况 下 根据 采用 视 在 光线 法 计算 得 到 的 
某 探测 方向 上 探测 器 接收 辐射 强度 分 布 , 采用 LSQR 
方法 在 介质 辐射 特性 参数 已 知情 况 下 重建 介质 三 维 
温度 分 布 。 并 在 此 基础 上 提出 LSQR-SPSO 混合 算 
法 ， 将 其 应 用 于 同时 重建 三 维 温度 场 及 辐射 特性 参 
数 。 结 果 表 明 ， 当 已 知 辐射 特性 参数 时 ， 无 论 有 无 
测量 误差 , LSQR 方法 都 可 较 好 地 重建 出 介质 的 三 允 
温度 场 ; 当 辐 射 特性 参数 未 知 时 ， 无 论 有 无 测量 误 
差 , 采用 LSQR-SPSO 混合 算法 都 可 较 好 地 同时 重建 
出 介质 的 吸收 、 散 射 系数 及 三 维 温度 场 ， 证 明 
LSQR-SPSO 混合 算法 具有 较 高 的 计算 精度 和 较 好 的 
稳定 性 ， 可 广泛 应 用 于 半 透 明 介质 辐射 特性 参数 和 
温度 场 重建 研究 。 另 外 ， 在 同时 重建 过 程 中 ， 介 质 
的 温度 场 比 辐射 特性 参数 更 容易 被 重建 得 到 。 
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